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The recent discovery of critical phenomena arising in gravitational collapse near the threshold of 
black hole formation is used to estimate the initial mass function of primordial black holes (PBHs) . 
It is argued that the scaling relation between black hole mass and initial perturbation found for a 
collapsing radiation fluid in an asymptotically flat space-time also applies to PBH formation in a 
Friedmann universe, indicating the possible formation of PBHs with masses much smaller than one 
horizon mass. Owing to the natural fine-tuning of initial conditions by the exponential decline of the 
probability distribution for primordial density fluctuations, sub-horizon mass PBHs are expected to 
form at all epochs. This result suggests that the constraints on the primordial fluctuation spectrum 
based on the abundance of PBHs at different mass scales may have to be revisited. 
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In this Letter, we consider the initial mass function 
(IMF) of primordial black holes (PBHs) formed in the 
process of gravitationally collapsing primordial density 
fluctuations in the radiation dominated phase of the early 
universe [|l],^ . Implications of the PBH number and mass 
spectrum with regard to their contribution to the cosmic 
density and the 7-ray background (via Hawking evapora- 
tion) have been employed to constrain the spectral index 
of the primordial fluctuation spectrum |^j4| . Two aspects 
play a central role in these studies: first, for each horizon- 
sized space-time region there exists a critical threshold 
value, 5c, for the density (or mass) contrast 6, separat- 
ing its further evolution between formation of a black 
hole {5 > 5c) and dispersion by pressure forces {5 < 5c) 
(we shall use the term "horizon" to denote the parti- 
cle horizon, rh ~ t). Comparing the Jeans and horizon 
lengths at the time when the collapsing region breaks 
away from Hubble expansion, one finds that 5c must be 
of order unity ||^ . The second key assumption relates to 
the final mass of the black hole, Mbh. It is commonly 
assumed that Mbh is approximately equal to the mass 
of the collapsing region and thus to the horizon mass 
at the epoch of formation, Afh. Nevertheless, detailed 
predictions for the PBH IMF have not previously been 
made. Based on a scaling relation discovered in gravita- 
tional collapse of various near-critical space-times, gener- 
alized to collapsing density perturbations in an Einstein- 
de Sitter universe, we are able to derive a universal, two- 
parameter PBH IMF, applicable when PBH number den- 
sities are dominated from fluctuations collapsing during 
one particular epoch. Here the two parameters in the 
PBH IMF carry all the information of the statistics of 
the initial density spectra and the perturbation shapes. 
We show that when the perturbation overdensity is suf- 
ficiently close to the critical overdensity for PBH forma- 
tion, 5c, the final mass of the resulting PBH may be an 



arbitrarily small fraction of the horizon mass, providing 
a conceptual difference to our current understanding of 
the process of PBH formation. 

It is possible that PBH formation is the only natural 
example for critical phenomena in gravitational collapse, 
a field of considerable interest in classical general relativ- 
ity that was previously believed to have no astrophysical 
application. Triggered by the intriguing results of Chop- 
tuik Q who demonstrated scaling and self-similarity in 
the gravitational collapse of a massless scalar field near 
the threshold of black hole formation, critical phenomena 
were studied for a number of different set-ups, includ- 
ing spherically symmetric radiation fluids , Yang-Mills 
fields 1^ , and axially symmetric collapsing gravity waves 
1^. In all cases, families of initial data quantified by a 
single generic parameter 5 were found to give rise to a 
scaling relation of the form 



Afbh(^) - if (<5-<5c)' 



(1) 



near the critical point for black hole formation, 5c. The 
specific choice of 5 is arbitrary since differentiable trans- 
formations of 5 leave (|l|) invariant, changing only the 
constant K to leading order Another noteworthy 

feature of near-critical solutions is the the appearance of 
discrete (scalar field collapse) or continuous (perfect fiuid 
collapse) self-similarities. 

Equation (0) is, in general, irrelevant for the forma- 
tion of astrophysical black holes. Degeneracy pressure of 
neutrons or electrons introduces intrinsic limiting mass 
scales of hydrostatic stability, such as the Chandrasekhar 
mass, violating the scale-free behavior indicated by equa- 
tion (|l|). Moreover, equation (^ is only valid in the 
immediate neighborhood of 5c, requiring a high degree 
of fine-tuning of the initial conditions which is unnat- 
ural under most circumstances. In PBH formation, on 
the other hand, it is expected that most regions collaps- 
ing to a black hole will have overdensities close to the 
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critical overdensity for PBH formation, 5c, owing to a 
steeply declining probability distribution for primordial 
density fluctuations. Typical cosmic initial conditions 
thus provide the fine-tuning of initial conditions, required 
for near-critical collapse. Further, with the exception of 
cosmological phase transitions that will not be considered 
here, the matter collapsing to PBHs is well described by a 
perfect fluid with radiation dominated equation of state, 
p = p/3, where p and p are pressure and energy density, 
respectively. Hence, the problem for PBH formation in 
radiation dominated cosmological epochs and the perfect 
fluid collapse studied numerically by Evans and Cole- 
man (EC) [0 differ only with regard to the background 
space-time. While canonical initial conditions for PBH 
formation involve curvature perturbations in an expand- 
ing Friedmann-Robertson- Walker (FRW) space-time, EC 
used initial data embedded in in an asymptotically sta- 
tionary space-time for their collapse simulations. 

In addition to their numerical simulations, EC found a 
self-similar solution to the equations of motion and grav- 
itation in the limit 6 6c- A self-similar ansatz reduces 
the spatial and temporal degrees of freedom to a single 
self-similar coordinate and thereby transforms the sys- 
tem of partial differential equations into ordinary ones. 
Demanding regularity at the center and along the ingo- 
ing acoustic characteristic, corresponding to the absence 
of a shock, the system of ODEs can be solved and the 
solution coincides well with their numerical results. As 
suggested by EC, the critical exponent of was sub- 
sequently derived by analyzing linear perturbations of 
the self-similar solution: Koike, Hara and Adachi pl| ] 
obtained 7 = 0.3558019 for a collapsing radiation fluid. 
Note that neither the self-similar solution nor the pertur- 
bation analysis rely on asymptotic flatness of the space- 
time; on the contrary, EC's self-similar solution is not 
asymptotically flat. As EC's solution converges neither 
to a flat stationary space-time nor to an exact FRW so- 
lution, it invariably breaks down at large radii for both 
asymptotic behaviors. 

The main reason to expect the emergence of self- 
similarity in near-critical gravitational collapse occuring 
in asymptotically FRW space-times is the separation of 
characteristic scales: Just as in the asymptotically sta- 
tionary case studied by EC, the solution forms an inter- 
mediate asymptotic between two widely separated length 
scales ||l2| . The scale rg of the fluid perturbation S at the 
onset of collapse is given by (5^^/-^rh 1^ if the initial per- 
turbation amplitude is evaluated at horizon crossing, tq 
can be identified with the transition from Hubble expan- 
sion of the asymptotic FRW space-time to the collapse- 
dominated region r < rg. On small scales, deviation 
from exact criticality leads to violation of self-similarity 
if r approaches ri ^ K\S — Sc\ The ratio ro/ri can 
be made arbitrarily large by chosing initial data close to 
the critical point. In the limit i5 ^ i5c, we therefore as- 
sume that gravitational collapse of a radiation fluid is 
well described by the self-similar solution of EC and 
the critical exponent 7 ~ 0.356 independent of the 



asymptotic behavior of the background space-time. We 
note that preliminary results of numerical simulations of 
the PBH formation process in the early universe confirm 
(1) the scaling relation and (2) the applicabilty of scal- 
ing for commonly assumed parameters of the statistics of 
pre-existing cosmic density fluctuations (see below) . The 
results of this numerical investigation will be presented 
elsewhere 10. 

Based on the arguments above we will henceforth em- 
ploy equation (|l|) for the masses of PBHs formed by col- 
lapsing primordial density perturbations slightly exceed- 
ing Sc, with an exponent 7 « 0.356 independent of initial 
perturbation shape. We assume a Gaussian probability 
distribution for density fluctuations entering the horizon, 

where a is the, possibly scale-dependent, root-mean- 
square fluctuation amplitude. Equation (H) allows us to 
compute the fraction of horizon-sized regions collapsing 
to PBHs at a given epoch jH) 

(3^ j\{S)dS^acxp(-^y (3) 

The upper integration limit reflects that US > 1, the 
collapsing space-time region corresponds to a separate 
closed universe instead of a black hole and the ap- 
proximation on the right hand side is valid to within a 
factor of a few for (t/6c < 0.2. It is noted that non- 
Gaussian effects may be important for Sc^ <J , but a 
Gaussian distribution suffices for the demonstration pur- 
pose of this work. 

In what follows we assume that cosmological PBH for- 
mation is dominated by perturbations of one particular 
length scale, defining a characteristic epoch of PBH for- 
mation by the time the perturbations cross into the hori- 
zon. Such an analysis should be adequate when cither 
the initial perturbation spectrum exhibits a peak on a 
given scale ]l6t , or PBH formation is most probable dur- 
ing a specific epoch by virtue of the equation of state 
. It is also approximately valid for blue initial pertur- 
bation spectra where PBH formation is most efficient on 
the smallest scale under consideration [T^ . 

With these assumptions, and approximating incor- 
rectly for the moment that the mass of the resulting PBH 
is Mh, we may compute the value of the PBH mass den- 
sity divided by the cosmic background density, 

fipbh.oid = (^) - M-W Mbh P{5) dS 
Po J 5, 

« /? for Mbh « Mh , (4) 

where the hat indicates that f2pbh is evaluated at the time 
of PBH formation. 

As a straightforward modification of flpbh 1 we use the 
continous distribution of PBH masses (|l[) in (|^) and re- 
evaluate the integral. Doing so, we implicitly assume that 
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is valid for S as large as unity; this need not neces- 
sarily be the case. However, the largest contribution to 
the integral comes from S ^ Sc owing to the exponential 
form of P{S), and thus our assumption is justified. The 
integrand rises steeply to a maximum at 

'^m ^l{Sc + ^47^' + 6^) - <5c + ^ + 0(^4) (5) 

close to the lower integration boundary, and the black 
hole mass at this point is 



Mbh(<5m) = K 



7(T 



2X7 



(6) 



with the dimensionless k defined by K = kM\,. The 
modified expression for f2pbh is thus 

J^pbh.ncw = Mi7' / Mu.{5) P[5) dS 

J 5^ 



^ka'+^-'cxpl^-^^ ^ka'''/3 , (7) 

where the integral was asymptotically expanded to first 
order. Equation shows that the average black hole 
mass produced at each epoch is approximately given by 
(m, since 



(Mbh) = p- 



Mbh((5) Pi6) d6 « fccr^'^Mh 



(8) 



We can now determine the PBH initial mass func- 
tion (IMF) when PBH number densities are dominated 
from formation during one particular epoch. The global 
PBH mass spectrum generally involves an integration 
over all epochs, a formidable problem, which will not 
be attempted here. We define the PBH IMF as the frac- 
tion dcj) of PBH number per logarithmic mass interval, 
normalized such that 



lnMbh(5=l) 



d(lnM^J 
This mass function is given by 



d(lnM^ 



bh^ 



(9) 
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d(lnMbh) 



r'pmihh)) 



dd 



d(lnAfbh) 



1 



2'Kl3a'-f 



"^bhexp 



'bh; 



2a2 



(10) 



where mbh is black hole mass in units of fcMh, and 
where we have used equation (|) for P{5). The PBH 
IMF of equation ( p^ ) has wider applicability than naively 
thought. Imagine PBH formation in the case of non- 
Gaussian statistics, in particular, when P{5) is differ- 
ent from equation (0). In this case one may search 



for a control parameter 5'{5) which renders P{5)d5/d5' 
Gaussian. Applying this transformation between con- 
trol parameters to equation ([|), one will obtain a form- 
invariant equation with modified constants K' and 
5'^, provided the limit of near-critical gravitational col- 
lapse still applies. Equation (|lo|) then defines a universal 
two-parameter family of PBH IMFs, applicable for many 
initial conditions, with the parameters K and 5c carrying 
all the information about the statistics of initial condi- 
tions and the shapes of perturbations. 

The mass function of equation (^0|) exhibits a maxi- 
mum at 



2 \ 7 



(11) 



which approximately equals the average black hole mass 
equation (H). PBHs in cosmological interesting numbers 
are formed during the evolution of the early universe for 
values of a/Sc ~ 0.1 — 0.2, provided Gaussian statistics 
holds Q]. Such values of a /5c yield typical volume col- 
lapse fractions in the range (3 k, 10~^ — 10~^^ and, de- 
pending on the epoch of formation, imply PBH number 
densities significantly contributing to the present mass 
density, or the 7-background. Inserting k « 3.3 found by 
EC |lgl, (j/5c w 0.15, and 5c ~ 1/3 [|, we find 



Mmax 
bh 



0.6Mh 



(12) 



It is not surprising that the maximum of the IMF at a 
fixed epoch coincides with the horizon mass to within an 
order of magnitude, since the latter determines the mass 
scale for collapse. However, depending on the value of 
a, a fraction of all PBHs formed at each epoch will have 
masses significantly smaller than Mb, implying a funda- 
mental conceptual difference between this work and pre- 
vious calculations. It was previously assumed that there 
exists a one-to-one correspondence between Mbh and red- 
shift z. Under this assumption, it was straightforward to 
relate Mbh to a single energy scale, i.e. microscopically 
small black holes only formed at very early times. Us- 
ing equation (|l|) instead, this simplification is no longer 
valid; the formation of black holes with a continuous IMF 
allows the formation of microscopic PBHs at all epochs. 
The formulation of observational constraints based on 
these results, such as constraints on the spectral index of 
initial density spectra, requires a detailed analysis of the 
PBH IMF integrated over all epochs which is beyond the 
scope of this Letter. 

The authors wish to thank S.E. Woosley for helpful 
discussions. 
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